	REAL(8) FUNCTION epsil(h0)
	implicit none
	real(8) h0, pi
      pi = 3.1415926
	   
	if (h0 >= 15./180.*pi.and.h0 < 20./180.*pi) epsil = 0.24
	if (h0 >= 20./180.*pi.and.h0 < 25./180.*pi) epsil = 0.22
	if (h0 >= 25./180.*pi.and.h0 < 30./180.*pi) epsil = 0.20
	if (h0 >= 30./180.*pi.and.h0 < 40./180.*pi) epsil = 0.18
	if (h0 >= 40./180.*pi.and.h0 < 50./180.*pi) epsil = 0.16
	if (h0 >= 50./180.*pi.and.h0 <= 60./180.*pi) epsil = 0.14
	if (h0 < 15./180.*pi) epsil = 0.30
		   
	RETURN
	END FUNCTION epsil     